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TITLE 

MEASURING METHOD FOR DECIDI NG DIRECTIO N TO A FLICKERING SOURC E 

— ■ 11 ' 

5 TECHNICAL FIELD 

The present invention relates to a method for deciding the 
direction to a flickering source in relation to a measurement 
point in an electrical network with alternating current with a 
network frequency with low-frequency amplitude variations from 
10 a flickering source. The present invention also relates to an 
arrangement t comprising means for carrying out the method. 

BACKGROUND ART 

During the production of electricity, generators are used that 
15 produce an alternating voltage around a certain frequency. The 
user of the produced electricity is aware of which frequency 
is involved and which voltage is provided on the electricity 
network. The user of the electricity network requires as clean 
electricity as possible, that is requires electricity that is 
20 well-defined at the specified frequency and the specified 
voltage. In normal electricity networks that supply towns and 
factories, it is, however, normal for the electricity network 
to comprise low- frequency voltage distortions. Voltage 
distortions arise when various devices are cyclically 
25 connected into the network, which devices comprise, for 
example, capacitances and inductances. 

Large power-hungry installations such as induction furnaces, 
compressors, lift motors, pumps, etc, help to increase the 
30 amount of voltage distortion in the electricity network. These 
loads often cause low-frequency (<25 Hz) and periodic 
fluctuations in the effective value of the voltage. The 
phenomenon is usually called flicker. The most noticeable 
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effect of flicker is the annoyance that can be experienced 
when the light intensity from incandescent lamps fluctuates in 
line with the voltage variations caused by the flicker. 
Studies have shown that people are particularly sensitive to 
5 light fluctuations with a repetition frequency in the range 
0.5 Hz to 25 Hz. At the level of maximum sensitivity (approx. 
9 Hz), the relative voltage change needs only to be 0.25% for 
people to be able to perceive that the light from the 
incandescent lamp flickers. Problems with flicker arise 
10 primarily in areas with very heavy industry (ironworks and 
paper-mills, etc) but can also arise in areas with weak 
electricity networks and in the vicinity of wind-power 
installations . 

15 In the event of a dirty electricity network comprising 
flicker, it is of interest to know where the flickering source 
is located. An electricity producer wants to show that it is a 
consumer that makes the electricity network dirty and can, in 
such a case, demand that the consumer pay a penalty charge or 

20 puts the problem right. A consumer, on the other hand, wants 
of course to show that it is not him causing the problem. In 
addition, the consumer wants to have value for money and thus 
wants to show that it is the electricity producer that is 
supplying dirty electricity. 

25 

Putting right the problem of flicker is often an expensive 
operation where all or parts of the network can need to be 
reconstructed in such a way that the internal impedance of the 
electricity network is reduced, for example by means of new 
30 and larger cables. Another way to put the problem right is to 
install counter measures against the flicker in question. Such 
counter measures are usually very expensive. An example of a 
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counter measure is "Static Var Compensation (SVC) " which 
dynamically controls changes in the system. 

Methods exist for determining the content of voltage flicker, 
5 which are described in the standard IEC-61000-4-15 . Measuring 
instruments that record voltage flicker according to this 
standard show the presence of voltage flicker by calculating 
and displaying the parameters Ifl, Pst and Pit. There is 
however nothing in the standard or in existing measuring 
10 instruments that shows the direction to the flickering source 
in relation to the measurement point. 

There is thus a great need for a method and an arrangement 
that can determine the direction of propagation of flicker. 
15 The direction of propagation of the flicker shows whether the 
flickering source is above or below a measurement point and 
can thus be of great use when the flickering source is to be 
traced. 

20 DISCLOSURE OF INVENTION 

The present invention intends to solve the problems that are 
described above of low-frequency flicker that gives rise to 
periodic fluctuations in the effective value of the voltage. 
In the following, such periodic fluctuations are called low- 

25 frequency variations or flicker. 

The problem is solved by means of a method for deciding the 
direction to a flickering source in relation to a measurement 
point in an electricity network with alternating current with 
30 a network frequency f c with low-frequency amplitude variations 
from the flickering source. The method is characterized in 
that an amplitude -modulated current signal and an amplitude- 
modulated voltage signal are recorded. Both the amplitude- 
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modulated current signal and the amplitude-modulated voltage 
signal are signal processed in such a way that the low- 
frequency amplitude variations in both the current and voltage 
are separated from the carrier wave in the form of a flicker 
5 component for current and a flicker component for voltage. The 
flicker component for current is thereafter multiplied by the 
flicker component for voltage in such a way that a product is 
produced. The product is processed in such a way that a 
flicker power n is obtained with a sign value that indicates 
10 in which direction a flickering source is located in relation 
to the measurement point . 

According to an embodiment of the invention, the sign value of 
the flicker power at the measurement point is negative when 
15 the flickering source is located below the measurement point 
and positive when the flickering source is located above the 
measurement point . 

The method has two important advantages: 
20 1. The flicker power at the measurement point is 

determined. By this means the flickering source can be traced. 

2. The method is also accurate when the measurement is 

carried out via voltage and current transformers installed in 

25 the network. The frequency spectrum of a signal with flicker 
consists of a carrier wave (for example 50 or 60 Hz) and of 
side-bands on both sides of the carrier wave, with frequency 
distance to the carrier wave corresponding to the flicker 
frequencies. The frequency spectrum of the whole signal 

30 package (carrier wave and the low- frequency signals) is thus 
in a narrow frequency band around the carrier wave, which 
means that measurement via existing current and voltage 
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transformers can be used as these are designed for highest 
accuracy around the network frequency. 

The invention can be carried out by means of a number of 
5 methods described below. 

In a first method (method 1) according to an embodiment of the 
invention, the method is characterized by the steps: 

recording of an amplitude-modulated current signal i (n) 
10 comprising signals that originate from the network frequency f c 
and the low-frequency amplitude variations in the current 
signal i (n) ; 

recording of an amplitude-modulated voltage signal u (n) 
comprising signals that originate from the network frequency f c 
15 and the low-frequency amplitude variations in the voltage 
signal u (n) ; 

- signal processing of the current signal i (n) in such a way 
that only the low-frequency amplitude variations remain in the 
form of a flicker component for the current signal i (n) ; 

20 - signal processing of the voltage signal u (n) such a way that 
only the low-frequency amplitude variations remain in the form 
of a flicker component for the voltage signal u(n); 

creation of a product by multiplication of the flicker 
component for current and the flicker component for voltage; 

25 - processing of the product in such a way that a flicker power 
FT is obtained with a sign value that indicates in which 
direction the flickering source is located in relation to the 
measurement point . 

30 An advantage of the first method is that it does not require 
extensive computer capacity, but can easily be implemented in 
a suitable device. 
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According to another embodiment of method 1, the method 
comprises : 

- the signal processing of the current signal i (n) comprising 
the steps : 

5 - creation of a first demodulated signal by demodulation 

of the current signal i (n) ; 

- filtering off of the signals that originate from the 
network frequency f c in the first demodulated signal in 
such a way that only the low- frequency variations remain 

10 in the form of the flicker component for current; 

- the signal processing of the voltage signal u (n) comprising 
the steps : 

- creation of a second demodulated signal by demodulation 
of the voltage signal ; 

15 - filtering off of the signals that originate from the 

network frequency f c in the second demodulated signal in 
such a way that only the low-frequency variations remain 
in the form of the flicker component for voltage. 

20 According to another embodiment of method 1, the method 
comprises the steps: 

- filtering off of the signals that originate from the network 
frequency f c in the first demodulated signal in such a way that 
only the low-frequency variations relating to the flicker 

25 component for current remain in the form of a flicker signal 
iLF(n) for current; 

- filtering off of the signals that originate from the network 
frequency in the second demodulated signal in such a way that 
the low-frequency variations relating to the flicker component 

30 for voltage remain in the form of a flicker signal U LF ( n) for 
voltage ; 
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- the product creating an instantaneous power signal n(n) by 
multiplication of the flicker signal I LF (n) for current and the 
flicker signal U L F(n) for voltage; 

- the product being processed by the creation of the average 
5 value of the instantaneous power signal Il(n) whereby the 

flicker power n is created with the sign value. 

According to yet another embodiment of method 1, the method 
comprises : 

10 - the first demodulated signal being created by square 
demodulation of the current signal; 

the second demodulated signal being created by square 
demodulation of the voltage signal . 

15 According to an embodiment of method 1, the filtering is 
carried out with a bandpass filter with a lower limit of 0.1 
Hz and an upper limit of 35 Hz. A preferred upper limit is, 
however, 2 5 Hz. 

20 In a second method (method 2) according to an embodiment of 
the invention, the method is characterized by the steps: 

recording of an amplitude-modulated current signal i (n) 
comprising signals that originate from the network frequency f c 
and the low-frequency amplitude variations in the current 

25 signal i (n) ; 

recording of an amplitude-modulated voltage signal u (n) 
comprising signals that originate from the network frequency f c 
and the low- frequency amplitude variations in the voltage 
signal u (n) ; 

30 - frequency analysis of the wave form of the voltage signal 
u(n) by means of an N-point DFT analysis (Discrete Fourier 
Transform) , which gives rise to a voltage vector U which 
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contains the frequency spectrum for the voltage signal u(n) in 
the form of N complex voltages; 

- frequency analysis of the wave form of the current signal 
i (n) by means of an N-point DFT analysis (Discrete Fourier 

5 Transform) , which gives rise to a current vector I which 
contains the frequency spectrum for the current signal i (n) in 
the form of N complex currents; 

- creation of a power vector P by means of the multiplication, 
element by element, of the voltage vector U and the current 

10 vector I ; 

- multiplication of the power vector P by a weighting vector W 
that eliminates the power component that originates from the 
network frequency, with the power vector P comprising partial 
powers P k concerning power components from the flickering 

15 source, 

- creation of a flicker power n with a sign value by means of 
summation of the partial powers P k , and 

- analysis of the sign value, with the sign value indicating 
in which direction from the measurement point the flickering 

20 source is to be found. 

In the second method, the voltage signal u (n) and the current 
signal i (n) are signal processed by means of the frequency 
analysis described above. The product given by method 1 

25 corresponds to the creation of the power vector P in method 2. 
The flicker components for current and voltage in method 1 do 
not have anything directly corresponding to them in method 2, 
but flicker components arise in the form of power components P k 
concerning power components from the flickering source after 

30 the power vector P has been multiplied, element by element, by 
the weighting vector W which comprises the elements w k . The 
processing described in method 1 corresponds to the summation 
of the partial powers P k in method 2. 
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According to an embodiment of method 2, the flicker power n is 
created by means of the following step: 

- summation of the partial powers P k by means of the formula: 

5 

An advantage of this method is that no demodulation residues 
arise . 

According to another embodiment of method 2, the flicker power 
10 n is created by means of the following steps 

- square demodulation x 2 of the voltage signal u(n); 

- square demodulation x 2 of the current signal i (n) ; 

calculation of the frequency spectrum of the square- 
demodulated voltage signal by means of an N-point DFT analysis 
15 (Discrete Fourier Transform) which gives rise to the voltage 
vector (U) ; 

calculation of the frequency spectrum of the square- 
demodulated current signal by means of an N-point DFT analysis 
(Discrete Fourier Transform) which gives rise to the current 
20 vector (I) ; 

- creation of the flicker power IT by means of summation of the 
partial powers P k which contribute to the flicker phenomenon by 
means of the formula: 

n =Z Re {^ £v>*v/;} 

25 where the elements wl k and w2 k replace W k and eliminate the 
power component that originates from the network frequency and 



9 



weight the correct amplitudes of the frequency components U k 
and I k/ in accordance with: 



for \<k<i 



0 



for k > i 



> 



0 



for k>i 



for \<k<i 
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where it is assumed that the low- frequency flickers are to be 
found in a frequency band up to and including tone i (0 < 



10 As shown above, both method 1 and method 2 give rise to the 
flicker power II according to the invention with a sign value 
that indicates in which direction a flickering source is 
located in relation to a measurement point. It is thus 
possible to signal process the current signal and voltage 

15 signal in both the time plane and the frequency plane in order 
to obtain the required flicker power EL 

The invention also relates to an arrangement comprising means 
for carrying out the methods described above . 



According to an embodiment of the invention, the arrangement 
comprises means for deciding the direction to a flickering 



f flicker ^ i) - 
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source in relation to a measurement point in an electricity- 
network with alternating current with a network frequency f c 
with low-frequency amplitude variations from the flickering 
source. The arrangement is characterized in that it also 
comprises : 

- means for recording an amplitude -modulated current signal 
i (n) comprising signals that originate from the network 
frequency f c and the low- frequency amplitude variations in the 
current signal i (n) ; 

- means for recording an amplitude-modulated voltage signal 
u(n) comprising signals that originate from the network 
frequency f c and the low- frequency amplitude variations in the 
voltage signal u(n); 

- means for signal processing the current signal i (n) in such 
a way that only the low-frequency amplitude variations remain 
in the form of a flicker component for the current signal 
i(n) ; 

- means for signal processing the voltage signal u(n) in such 
a way that only the low-frequency amplitude variations remain 
in the form of a flicker component for the voltage signal 
u (n) ; 

means for creating a product by multiplication of the 
flicker component for current and the flicker component for 
voltage; 

- means for processing the product in such a way that a 
flicker power IT is obtained with a sign value that indicates 
in which direction the flickering source is located in 
relation to the measurement point. 

According to another embodiment of the invention: 

- the means for signal processing of the current signal (i(n) 
comprises : 
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- means for creating a first demodulated signal by means 
of demodulation of the current signal i (n) ; 

- means for filtering off the signals that originate from 
the network frequency f c in the first demodulated signal 
in such a way that only the low- frequency variations 
remain in the form of the flicker component for current; 

- the means for signal processing of the voltage signal u(n) 
comprises : 

- means for creating a second demodulated signal by means 
of demodulation of the voltage signal; 

- means for filtering off the signals that originate from 
the network frequency in the second demodulated signal in 
such a way that only the low-frequency variations remain 
in the form of the flicker component for voltage. 

According to another embodiment of the invention, the 
arrangement comprises : 

- means for frequency analysis of the wave form of the voltage 
signal u (n) by means of an N-point DFT analysis (Discrete 
Fourier Transform) , which gives rise to a voltage vector U 
which contains the frequency spectrum for the voltage signal 
u(n) in the form of N complex voltages; 

- means for frequency analysis of the wave form of the current 
signal i (n) by means of an N-point DFT analysis (Discrete 
Fourier Transform) , which gives rise to a current vector I 
which contains the frequency spectrum for the current signal 
i (n) in the form of N complex currents; 

means for the creation of a power vector P by the 
multiplication, element by element, of the voltage vector U 
and the current vector I; 

- means for the multiplication of the power vector P by a 
weighting vector W that eliminates the power component that 
originates from the network frequency, with the power vector P 
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comprising partial powers P k concerning power components from 
the flickering source, 

- means for the creation of a flicker power n with a sign 
value, by summation of the partial powers P k , and 
5 - means for analysis of the sign value, with the sign value 
indicating in which direction from the measurement point the 
flickering source is to be found. 

BRIEF DESCRIPTION OF DRAWINGS 
10 The invention will be described below in the form of a number 
of embodiments with reference to a number of figures, in 
which: 

Figure 1 shows an equivalent two-terminal network for an 
15 electricity network according to the invention; 

Figure 2a shows variations in the effective value U L of the 
voltage; 

20 Figure 2b shows variations in the effective value I L of the 
current ; 

Figure 3 shows a frequency spectrum for an amplitude-modulated 
voltage signal with only one low-frequency component; 

25 

Figure 4 shows a signal flow chart for a measurement method 
according to an embodiment of the invention; 

Figure 5 shows a signal flow chart for a measurement method 
30 according to a second embodiment of the invention; 

Figure 6 shows a signal flow chart for a measurement method 
according to yet another embodiment of the invention; 
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Figure 7 shows an amplitude characteristic for a bandpass 
filter that is used in the embodiment described in association 
with Figure 6; 

5 

Figure 8 shows an outline block diagram of an arrangement that 
can be used for the methods described in the different 
embodiments ; 

10 Figure 9 shows schematically a network comprising a flickering 
source Fl, a load LI, a generator G for generating an 
alternating voltage, and 

Figure 10 shows schematically a diagram of the flicker power II 
15 for a number of sampling points n. 

PREFERRED EMBODIMENTS 

Both the general theory and a number of embodiments will be 
described below. The general theory is necessary for an 
20 understanding of the embodiments described below. In the 
equations that are described, a point between two letters in 
an equation indicates a multiplication of vectors, element by 
element . 

25 As mentioned above, the invention relates to a method for 
deciding in which direction a flickering source is located in 
relation to a measurement point. To determine the direction of 
low-frequency variations such as flicker, the wave form of the 
voltage and current must be recorded in the phase or phases 

30 where the direction is to be determined. Thereafter, the 
recorded information is signal processed according to any one 
of the embodiments that are described in the following signal 
flow charts. The result of the signal processing is a flicker 
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power n with a sign value. The sign value indicates in which 
direction the flickering source is located in relation to a 
measurement point . 



5 Figure 1 shows an equivalent two- terminal network for an 
electricity network according to previously known technology. 
The electricity network can be divided schematically into 
three parts which are commonly called generator 1, 
transmission line 2 and load 3. When a load 3 is connected, a 

10 current I will flow in the circuit. This leads to a drop in 
voltage U ZT arising across the internal impedance Z T so that 
the voltage U L across the load drops. If the load is connected 
in and disconnected periodically (cyclically) , the current I 
will also increase and decrease cyclically, so that the 

15 voltage U L decreases and increases cyclically (provided that 
the generator voltage U G remains constant) . 

Figure 1 shows a measurement point 17 and a point 18 that 
marks a point above, that is upstream of, the measurement 
20 point 17, and a point 19 that marks a point below, that is 
downstream of, the measurement point 17 . The terms above and 
below the measurement point 17 are important, as it is 
necessary to indicate in which direction the flickering source 
is located in relation to the measurement point. 

25 

According to the invention, a recording is taken at the 
measurement point of an amplitude-modulated current signal 
i (n) comprising signals that originate from a network 
frequency f c and low- frequency amplitude variations in the 
30 current signal i (n) . The low- frequency amplitude variations 
originate from the flickering source. In addition, a recording 
is carried out of an amplitude-modulated voltage signal u(n) 
comprising signals that originate from a network frequency f c 
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and the low-frequency amplitude variations in the voltage 
signal u(n). Also here, the low-frequency amplitude variations 
originate from the flickering source. 

5 Figures 2a and 2b show outline diagrams of variations in the 
effective value Urm S of the voltage for the voltage U L and 
effective value I 

rms of the current for the current voltage I . 
According to the invention, the changes in the effective 
values for current and voltage Urms and Irms relate to the 

10 instantaneous changes in the voltage and current that arise 
due to the low-frequency amplitude variations from the 
flickering source. Figures 2a and 2b do not show signals that 
originate from the network frequency. It is, however, known 
that the low-frequency amplitude variations modulate the 

15 network frequency f c . Any changes that take place in the 
current and voltage and that originate from the network 
frequency are negligible. The effective value changes Urm S and 
Irms thus reflect the instantaneous changes in current and 
voltage that can be derived from the low-frequency variations. 

20 In Figures 2a and 2b, from the point of view of the invention, 
it is thus only of interest to show the variations in the 
effective value Urms of the voltage and the effective value Irm S 
of the current. The variations in the effective value Urms of 
the voltage and the effective value Irms of the current depend, 

25 according to the method according to the invention described 
below, on whether the flickering source is located above or 
below the measurement point 17 . 

Figure 2a shows the situation when the flicker power is 
30 propagated from the load towards the generator, with the 
changes in the effective values for current Irms and voltage 
Urms taking place instantaneously and in antiphase. When the 
load increases, the current I increases, whereby the drop in 
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voltage across U ZT increases, which means that U L reduces 
instantaneously . 

Figure 2b shows the situation when the flicker power 
5 propagates from the generator towards the load, with the 
changes in the effective values for current Irms and voltage 
Urms taking place instantaneously and in phase. An increase of 
U G gives an instantaneous increase of the current I, which 
gives an instantaneous increase of U L , which thus gives a 
10 simultaneous change of current and voltage. 

In a theoretical consideration of flicker, it is expedient to 
let the current signal and voltage signal u(t), i(t) be 
described as an amplitude modulation. The mathematical 
15 expression for such signals is given by [1] and [2] , in 
accordance with: 



The signals consist partly of a carrier wave U c cos(co c t) and 
I c cos(co c t) respectively, and the low- frequency flicker signals 
20 U mk cos(co k t) and I mk cos (co k t ) respectively. The index c relates to 
the carrier wave's component and the index m relates to the 
low-frequency variation's component. The index k relates to 
the index for a summation across values k=l to N. 

25 Figure 3 shows a frequency spectrum for an amplitude-modulated 
voltage signal with only one low- frequency component (single- 
tone modulat ion) . The low- frequency component relates here to 
a low- frequency tone that gives rise to the variations in 




[1] 




[2] 
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amplitude of both current and voltage. By low-frequency tone 
is meant here a low- frequency signal. 

Figure 3 shows the frequency bands to which the low- frequency 
variation gives rise. In Figure 3, f c relates to the carrier 
frequency and f m to the frequency of the low- frequency 
variation. The frequency components with single- tone 
modulation are to be found at the frequencies f c , f c -fm and 

f c + f m • 

Figure 3 also shows the amplitude U c for the voltage of the 
carrier frequency f c and the amplitude U m /2 for the voltage 
component from the low-frequency variation f c . 

15 Figure 3 shows that a frequency spectrum of the modulating 
tone forms an upper side-band f c +fm and a lower side-band f c -fm 
with half the original amplitude U m in each side-band and 
located at a frequency distance from the carrier frequency f c 
corresponding to the modulating frequency f m . Low frequencies 

20 such as the low-frequency flicker give rise to side-bands 
close to the carrier frequency. The higher the frequency of 
the flicker, the larger the frequency distance of the 
modulating frequency from the carrier frequency. 

25 The flicker power, II, is the power that has its origin in the 
modulating tones and can be expressed as: 

n = i jf Z U «* cos K' + A )TZ^* c°s(<V + ) \t = { orthogonally } = 

1 oU=i A*=» J £ 3 j 
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The formula [3] shows that, after multiplication and 
integration, the individual low- frequency tones in current and 
voltage create the flicker power. In order to determine this 
power, the amplitude and phase of the low-frequency signals in 
5 voltage and current must be known and the low- frequency 
signals must be able to be extracted from the signal package 
[1] and [2] . This can be carried out in several different ways 
that lead to the same result. Examples of different methods 
are described below as different embodiments of the invention. 

10 

Two embodiments will be described below that are based on the 
frequency spectrum of the sampled wave forms u [n] and i [n] 
being determined by carrying out an N-point DFT analysis (DFT 
= Discrete Fourier Transform) . In the practical case, an FFT 

15 analysis (FFT = Fast Fourier Transform) can be used, which is 
a calculation algorithm that provides the same information as 
the Fourier Transform. The output data from the analysis is 
two complex-value voltage and current vectors U and I, which 
contain frequency spectrums for u [n] and i [n] in the form of k 

20 = N complex voltages U k and currents I k . The frequency 
resolution that is obtained depends upon the selected sampling 
frequency (f s ) and the number of samples N that are included in 
the calculation, in accordance with the relationship Af = f g 
/N. 

25 

Mathematically, the vectors U and I are expressed: 

U = [c/ 1 ,C/ 2J t/ 35 ...(/„] = \u i \^P l p^/3 2 p,\Z0 i ,...p N \Z/} N J [4] 

I = \h J 2 J 3 »•• • J N J = Iki > | 7 2 K«2 > \h \^<* 3 ,...,\I N \Za N J [5] 
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The complex power, S, is created on the basis of the vectors U 
and I and contains the power quantities P (average power) and 
Q (reactive power) per tone in accordance with: 

[6] 

) 

5 The active power Pi is obtained as the real part of Si, in 
accordance with: 

P,=Re{S l }=Re{Pt+jQ l } [7] 

The vectors U, I and S contain the frequency spectrum from the 
angular frequency f = 0 Hz to f = f s Hz. There is often a 
10 desire to weight the information and/or zero certain 
frequencies. A simple and effective way of doing this is to 
introduce the weighting vector W, in accordance with: 

w =[wi,w 2 ,w 3 ,...,w w J 

[8] 

15 

The weighing factor W contains elements, w k , that contain 
constants that are multiplied, element by element, by U, I and 
S. If the voltage vector U is to be zeroed for frequencies 
above half the sampling frequency, the following operation is 
20 carried out (the point notation in [9] relates to 
multiplication element by element) : 
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where 



[9] 



N ' 

1 for l<k< — 
2 



0 for —<k<N 
2 



U m od is the modified voltage vector containing the frequency 
spectrum up to half the sampling frequency. Above this, the 
elements of the voltage vector are zero. 

The weighting vector's elements, w k , can be selected so that 
the required result is obtained. For example, a filter 
characteristic can be achieved by selecting suitable values of 
the elements w k . In addition, different weighting vectors can 
be introduced for U, I and S with the aim of achieving the 
required result. 

The flicker power, II, is obtained by summation of the powers, 
P k , which contribute to the flicker phenomenon. This is to say: 



Figure 4 shows a signal flow chart for a measurement method 
according to an embodiment of the invention for a first 
measurement method. The calculation of the flicker power n is 




[10] 
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carried out by frequency analysing the wave form of the 
voltage and current. 

The input signals for the signal flow chart in Figure 4 
consist of the two input signal vectors, u [n] and i [n] , which 
contain the sampled wave forms for voltage and current. The 
distance in time between two indices (n and n+1) in the input 
signal vectors corresponds to l/f s/ where f 8 relates to the 
sampling frequency. 

By utilizing the previously described theory for the signal 
flow chart that is described in Figure 4, the flicker power 
can be calculated. Figure 4 shows two blocks designated DFT 
that refer to a transformation of the current signal and 
voltage signal, i (n) and u(n) respectively, by means of a so- 
called Discrete Fourier Transform. The signal that comes out 
from the DFT for the voltage signal u(n) is designated U and 
comprises the components U k Z(3 k which relate to a number of 
values for different angles p k/ where the index k indicates 
which point is referred to in the DFT. The signal that comes 
out from the DFT for the current signal i (n) is designated I 
and comprises the components I k Zcp k which relate to a number of 
values for different angles cp k , where the index k indicates 
which point is referred to in the DFT. Figure 4 shows a 
multiplier 3 in which the vectors U, I and W are multiplied. 
The output signal from the multiplier is designated P and 
relates to a power vector comprising components partial powers 
P k that describe the active power in accordance with [7] as the 
real part of the complex power S in accordance with [6] and 
[7] . 

P k = Re{0.5-w k - 1 U k I Zp k . | l k | Z-(p k } 
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[10a] 

Figure 4 shows, in addition, a summation point Z that 
5 indicates the creation of the flicker power Ft by summation of 
the active powers P k in accordance with [10a], where: 

n=£ Re{0.5w k -|u k |z3k-llklz-<p k } 

[10b] 

10 The frequency spectrums of the wave forms u [n] and i [n] are 
calculated by an N-point DFT analysis. Thereafter, the flicker 
power II is calculated by utilizing [10] . The content of the 
weighting vector W is decisive for the final result. Elements 
in W are to be selected so that the power in the carrier wave 

15 is zeroed and so that only the tones that originate from 
flicker are included. This can be carried out in three ways. 
The information is used from either 1) the upper side-band, 2) 
the lower side -band or 3) there is a summation per frequency 
of the information in both side-bands, and the power is 

20 thereafter calculated by utilizing [10] . 

In order to exemplify how the calculations are carried out, we 
assume that the carrier wave signal is found in element M of 
the vectors U and I. We also assume that the low-frequency 
25 flicker tones in the vectors U and I are found in the elements 
k with the index: 

M -i<k <M and M <k<M + i 
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The number of flicker tones is given by the constant i. If the 
choice is made to utilize the information in the lower side- 
band for the calculation of the flicker power, the elements w k 
in the weighting vector W are selected, in accordance with 



0 for \<k<M-i 
2 for M-t<k<M 
0 for M <k<N 



The flicker power is thereafter created by utilizing the 
formula [10] . 

10 If the choice is made to utilize the information in the upper 
side-band for the calculation of the flicker power, the 
elements w k in the weighting vector W are to be selected, in 
accordance with 



0 for l<k<M 
2 for M<k<M + i 
0 for M+i<k<N 
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The flicker power is thereafter created by utilizing the 
formula [10] . 

If the choice is made to utilize the information in both side- 
bands for the calculation of the flicker power, the elements w k 
in the weighting vector W are selected, in accordance with 



0 for \<k<M-i 

1 for M-i<k<M and M<k<M+i 
0 for k = M 



The flicker power n is thereafter created by means of the 
formula 
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[10c] 



Figure 5 shows a signal flow chart for a second measurement 
method according to another embodiment of the invention, where 
the calculation of the flicker power is carried out by means 
of square demodulation of voltage and current, after which a 
DFT analysis of the low- frequency signals is carried out. 

The method is reminiscent of the first measurement method 
described in association with Figure 4, with the difference 
that, before the frequency spectrum is calculated, the signals 
are square demodulated, which is designated X 2 in Figure 5. By 
this means, the low- frequency signals are separated from the 
carrier wave, in accordance with the expression below: 



Square demodulation means that two signal packages are 
created, separated in frequency. One signal package consists 
of a direct component, the modulating frequencies and the 
modulating signals squared. The latter is an unwanted 






[11] 
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combination product. The second signal package contains the 
same terms but centred in frequency around double the carrier 
wave frequency. For calculation of the flicker power, only the 
terms are included that are marked by double underlining in 
5 [11] . One term contains the low-frequency flicker tones 
multiplied by the carrier wave amplitude and the other term is 
the modulated signal squared. In terms of size, the first term 
is much larger than the second, which means that the squared 
term affects the result very little. Corresponding expressions 
10 are obtained when squaring the current signal i (t) . 

The input signals for the signal flow chart in Figure 5 
consist of the two input signal vectors, u [n] and i [n] , which 
contain the sampled wave forms for voltage and current. The 

15 distance in time between two indices (n and n+1) in the input 
signal vectors corresponds to l/f s . In X 2 each element in the 
input signal vectors is squared and constitutes input data for 
the N-point Discrete Fourier Transform, DFT . The output signal 
from the respective DFT is complex-value vectors U and I 

20 containing frequency spectrums for u [n] and i [n] (both 
amplitude information and phase information) with a frequency 
resolution of Af = f s /N. The content in U and I is expressed 
by: 

U =[^ 1 ,L/ 2J ^,...^J=[Lr i |ZA,|C/ 2 |Z^ 2 ,|f/ 3 |ZA,..,|^|^,] 

I = [l l ,I 2 ,I 3 ,...,I Nt ]= ll l \Za l ,\l 2 \Za 2 ,\l 3 \Za 3 ,...,\l„\Za N ,] 
25 The complex power, S, is calculated by means of: 
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s=^w\ k -u wi k r 



[12] 



5 The sought flicker power n is then obtained by means of the 
formula 



[13] 



The elements wl k and w2 k zero the frequency components that do 
not cause flicker and weight the correct amplitudes of the 
frequency components U k and I k , in accordance with 



1 



wl k =\ U, 



— for \<k<i 
0 for k>i 
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[14] 



w2 A = 



— for \<k<i 

^ c 

0 for k > i 
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[15] 

In [13] , [14] and [15] , it is assumed that flicker tones are 
found in a frequency band up to and including tone i (0 < 
f flicker < i) . The correct value for the weighting factors can be 
identified by studying [11] . The method is not limited to the 
weighting factors that are given above, other constants in the 
weighting vector being able to be used to obtain the required 
filter effect. 

Figure 6 shows a signal flow chart according to yet another 
embodiment of the invention for a third measurement method, 
where the flicker power is calculated by square demodulation 
of voltage u(n) and current i (n) where the low- frequency tones 
are filtered by the use of bandpass filters 2A and 2B 
respectively. 

Instead of using DFT analysis and weighting vectors as 
described in the two previous embodiments described in 
20 association with Figures 4 and 5, in this measurement method 
bandpass filters 2A, 2B are used for filtering the low- 
frequency flicker tones. 

The input signal vectors u [n] and i [n] contain the sampled 
25 wave forms for voltage and current. In the components 1A and 
IB, the carrier wave (signal from the network frequency) and 
the low-frequency flicker tones are separated, similarly to 
the measurement method described in association with Figure 5, 
by each sample being squared. Only the low- frequency flicker 
30 tones are allowed to pass through the bandpass filters 2A and 
2B. As input signals to the multiplier 4, there are thus only 
the low- frequency flicker tones in voltage u LF [n] and current 
i L F[n] . The output signal from the multiplier 4 is the 
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instantaneous flicker power Ft [n] = u LF [n] x i L F[n], The flicker 
power IT is obtained by the integrator 5 creating an average 
value for the instantaneous power IT(n) . This can be carried 
out with a digital filter in the form of, for example, a low- 
5 pass filter. 

The bandpass filters 2A and 2B are dimensioned to obtain a 
lower limit frequency of 0.1 Hz and an upper limit frequency 
of 25 Hz. Alternatively, the bandpass filters can be used that 
10 are defined in the standard IEC 61000-4-15 that describes a 
flicker algorithm . 

The transmission function for the filter in the IEC standard 
is expressed as 



kco x s 



s 2 +2As + co 2 



(s + l/a> 3 )(s + l/& 4 ) 



15 



[16] 

The coefficients in [16] are to have values in accordance with 
20 the following table. 



k = 1.74802 


X = 27C-4. 05981 


©1=271-9. 15494 


C0 3 = 27l-1 .22535 


C0 2 = 27l-2 . 27979 


C0 4 = 27l-21 . 9 



Table 1. 



25 Figure 7 shows the amplitude characteristic for the bandpass 
filter with transmission function in accordance with [16] and 
coefficients in accordance with Table 1. 
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It is also possible to select bandpass filters with a 
different characteristic to that described in [16] . For 
example, an Mth order Butterworth or Chebishev filter can be 
5 selected. 

Figure 8 shows a block diagram of the hardware for the 
instrument. The instrument is constructed around a signal 
processor 7 which administers the measurement information and 
carries out necessary calculations according to the 
measurement methods that are described in the Figures 4-6 and 
9. The signal processor 7 also controls the sampling process 
for the conversion of analogue signals to digital signals (the 
A/D conversion) . The recorded signals, that is the wave form 
of current and voltage, are obtained either from current 
transformers and voltage transformers in the network or from 
measurement sensors that are to be found in the instrument. 
The signal processor can be a computer or a logic circuit, or 
some other suitable device for controlling devices and for 
signal processing of signals. 

Figure 8 shows a signal -conditioning device 8 for the recorded 
voltage signal. The measured voltage is signal -matched by 
means of a resistive voltage division, whereby the correct 
25 input signal level is obtained for the following step which is 
an anti-alias filter 10. 

Figure 8 also shows a signal - conditioning device 9 for the 
recorded current signal. The signal level for the current 
30 channel is matched to the instrument, either via a low-ohm 
shunt, with the voltage drop across the shunt being amplified 
and becoming the input signal level for the following step 
which is an anti-alias filter 11, or alternatively, the 
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current signals can be obtained from current tongs connected 
to the instrument . 

The task of the anti-alias filters 10, 11 is to prevent 
5 weighting distortion which arises if the recorded signal has a 
frequency content that exceeds half the sampling frequency 
(see theory for the sampling theorem) . The anti-alias filters 
can be implemented in the form of an analogue Sallen-Key low- 
pass filter according to Figure 7 and have a limit frequency 
10 corresponding to half the sampling frequency. 

The level -matched and filtered signals are sampled in sampling 
devices 12, 13 with a sampling frequency, for example 64 00 Hz. 
The digital raw data in the form of sampled amplitude values 
15 is saved in a measurement memory 14 and later constitutes the 
input data for the measurement methods described above . 

Software that controls the signal processor 7 is found in a 
program memory 15. The finished result, that is the flicker 
20 power with sign value and internal impedance, can be shown 
both numerically and graphically in a display device 16. The 
display device can be any known device for numerical and 
graphical display, for example a VDU. 

25 Figure 9 shows schematically a network comprising a flickering 
source Fl, a load LI, a generator G for generating an 
alternating voltage. In Figure 9, the current direction for 
the current I is shown by an arrow in the connecting lines. 
The current goes from the generator G to the flickering source 

30 Fl and the load LI. Figure 9 also shows a first measurement 
point Ml. In association with Ml, there is marked a point 
above 18 the measurement point Ml and a point below 19 the 
measurement point Ml. Figure 9 also shows a second measurement 
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point M2 . In association with M2 there is marked a point above 
18 the measurement point M2 and a point below 19 the 
measurement point M2 . The flickering source Fl emits into the 
network a low-frequency amplitude variation that propagates in 
5 the direction shown by a solid-line arrow 20 in Figure 9. 

In the method according to the invention described above, the 
flicker power II indicates in which direction the flickering 
source is located in relation to a measurement point. 
10 According to an embodiment of the invention, the sign value of 
the flicker power is negative when the flickering source is 
located below 19 the measurement point and positive when the 
flickering source is located above 18 the measurement point. 

At the first measurement point Ml, a negative sign is thus 
obtained for the flicker power II, as the flickering source Fl 
is located below 19 the first measurement point. This is 
because the low- frequency variations in current and voltage 
are in antiphase at the first measurement point Ml. 

At the second measurement point M2 , on the other hand, a 
positive sign is obtained for the flicker power n, as the 
flickering source Fl is located above 18 the second 
measurement point M2 . This is because the low- frequency 
variations in current and voltage are in phase at the second 
measurement point M2 . 

Figure 10 shows schematically a diagram of the flicker power II 
for a number of sampling points n during a certain period of 
30 time when recordings were made of the modulated current 
signals and voltage signals i (n) , u(n) . 
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Figure 10 shows a first curve Kl on the negative lower part of 
the diagram. The diagram also shows a second curve K2 on the 
upper part of the diagram. The first curve Kl corresponds to a 
power signal which, after the creation of the average value, 
5 gives rise to a flicker power n with a negative value and thus 
corresponds to the flicker power n that was obtained at the 
first measurement point Ml in Figure 9. The second curve K2 
corresponds to a power signal which, after the creation of the 
average value, gives rise to a flicker power IT with a positive 
10 value and thus corresponds to the flicker power II which was 
obtained at the second measurement point M2 in Figure 9. 

Figure 10 also shows that Kl inverted about the sampling axis 
n corresponds to K2 . This has been demonstrated in experiments 

15 and the explanation is that the flicker power IT changes sign 
when the propagation of the low- frequency flicker is changed 
from going against the power direction of the basic tone to 
going with the power direction of the basic tone, or vice 
versa. It is thus necessary to know in which direction the 

20 generator and load respectively are to be found in relation to 
the measurement point in order for it to be possible to 
interpret the sign value in the correct way. The part of the 
diagram that shows the period between zero and where Kl and K2 
respectively start, shows a period of time when the flickering 

25 source Fl is not connected in. 
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